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ABSTRACT 


The Sloan Lens ACS (SLACS) Survey is an efficient Hubble Space Telescope (HST ) Snapshot imaging survey for 
new galaxy-scale strong gravitational lenses. The targeted lens candidates are selected spectroscopically from the 
Sloan Digital Sky Survey (SDSS) database of galaxy spectra for having multiple nebular emission lines at a redshift 
significantly higher than that of the SDSS target galaxy. The SLACS survey is optimized to detect bright early-type 
lens galaxies with faint lensed sources in order to increase the sample of known gravitational lenses suitable for 
detailed lensing, photometric, and dynamical modeling. In this paper, the first in a series on the current results of our 
HST Cycle 13 imaging survey, we present a catalog of 19 newly discovered gravitational lenses, along with nine 
other observed candidate systems that are either possible lenses, nonlenses, or nondetections. The survey efficiency 
is thus >68%. We also present Gemini 8 m and Magellan 6.5 m integral-field spectroscopic data for nine of the 
SLACS targets, which further support the lensing interpretation. A new method for the effective subtraction of 
foreground galaxy images to reveal faint background features is presented. We show that the SLACS lens galaxies 
have colors and ellipticities typical of the spectroscopic parent sample from which they are drawn (SDSS luminous 
red galaxies and quiescent MAIN sample galaxies), but are somewhat brighter and more centrally concentrated. 
Several explanations for the latter bias are suggested. The SLACS survey provides the first statistically significant 
and homogeneously selected sample of bright early-type lens galaxies, furnishing a powerful probe of the struc- 
ture of early-type galaxies within the half-light radius. The high confirmation rate of lenses in the SLACS survey 
suggests consideration of spectroscopic lens discovery as an explicit science goal of future spectroscopic galaxy 


surveys. 


Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation — 
galaxies: structure — gravitational lensing 


1. INTRODUCTION 


In the currently favored cosmological scenario, the matter 
content of the universe is dominated by a cold and dark compo- 
nent of unknown particle species whose only significant interac- 
tion with the smaller baryonic matter fraction occurs through the 
gravitational force. This cold dark matter (CDM) picture is most 
strongly required by observations on the largest scales (Spergel 
et al. 2003; Tegmark et al. 2004; Percival et al. 2002). The CDM 
scenario holds that galaxies form within the potential wells of ex- 
tended dark matter halos that began their collapse while baryonic 
matter was still ionized (White & Rees 1978; Blumenthal et al. 
1984). This view is supported by direct evidence for dark matter 
in disk galaxies from the observation of rotational velocities that 
remain approximately constant out to radii at which the stellar 
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galactic component makes a diminishing contribution (e.g., Rubin 
et al. 1980; van Albada & Sancisi 1986). 

Unlike disk galaxies, early-type galaxies (E and SO) are gen- 
erally pressure supported and lack bright kinematic tracers at 
large radii such as H 1. As such, their kinematics are more difficult 
to measure and interpret in terms of mass density profiles. The 
density structure of early-type galaxies is nevertheless of great 
interest for numerous reasons. First, their structure is the fossil 
record of their formation and evolution processes (e.g., Wechsler 
et al. 2002; Zhao et al. 2003; Loeb & Peebles 2003; Gao et al. 
2004). Hierarchical CDM galaxy formation theories hold that 
early-type galaxies are built through the merging of late types 
(Kauffmann et al. 1993; Baugh et al. 1996), which should have 
predictable consequences for the structure of the merger prod- 
ucts. The most stringent test of these theories will require pre- 
cise observational measurements of early-type mass profiles. For 
example, detailed measurement of the structure of high surface 
brightness early-type galaxies will enable quantitative tests of 
the CDM theory on scales on which baryonic and radiative pro- 
cesses have significant effects on the structure of the host dark 
matter halo (e.g., through adiabatic contraction, e.g., Flores et al. 
1993; Blumenthal et al. 1986; Mo et al. 1998; Gnedin et al. 
2004), altering it significantly relative to the form expected to 
result from pure collisionless dark matter collapse (e.g., Navarro 
etal. 1996; Moore et al. 1998). Second, early-type galaxies exhibit 
great regularity in their photometric, spectroscopic, and kinematic 
properties, as described, e.g., by the well-known fundamental 
plane (FP) relation between velocity dispersion, effective radius, 
and surface brightness ( Djorgovski & Davis 1987; Dressler et al. 
1987). The tilt of the FP relative to the simple expectation of the 
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virial theorem can be understood in terms of a dependence of gal- 
axy structure and mass-to-light ratio on total galaxy mass (Bender 
etal. 1992; Ciotti et al. 1996; Bertin et al. 2002; Trujillo et al. 2004). 
However, further constraints on the mass structure of early-type 
galaxies are needed in order to distinguish between the various 
effects of differing stellar populations (the explanation given by, 
e.g., Gerhard et al. 2001), differing density profiles (structural 
nonhomology or weak homology), and differing dark matter frac- 
tions (the explanation of, e.g., Padmanabhan et al. 2004) in giving 
rise to the FP. 

Despite great progress, observational results remain uncer- 
tain because of the small number of individual galaxies suitable 
for study. Stellar dynamical measurements of local early-type 
galaxies (e.g., Bertin et al. 1994; Franx et al. 1994; Gerhard et al. 
2001; Cappellari et al. 2006), the statistics of early-type gravi- 
tational lens galaxies (e.g., Kochanek 1996; Chae 2003; Rusin 
et al. 2003; Rusin & Kochanek 2005), and combined lensing and 
dynamical measurements of the few systems amenable to such 
study (Koopmans & Treu 2002, 2003; Treu & Koopmans 2002, 
2003, 2004; hereafter collectively KT) generally argue for the 
presence of a significant amount of dark matter, even on the scale 
of the half-light radius, leading to an approximately linear in- 
crease of enclosed mass with radius and thus to a flat equivalent 
rotation curve as in disk galaxies. However, in some cases the 
presence of dark matter is not required, and the observed kine- 
matics can be reproduced with a constant mass-to-light ratio 
(e.g., Bertin et al. 1994; Romanowsky et al. 2003), correspond- 
ing to a declining rotation curve. Furthermore, Kochanek (2003) 
has pointed out an apparent conflict between isothermal mass 
profiles (i.e., flat rotation curves) and some gravitational lens 
time delays under the assumption of Hy ~ 70 km s7! Mpc™!. 
Due to this persistent uncertainty about the structure and diver- 
sity of early-type galaxies (see also Kochanek 2005), it is im- 
portant to collect precision data for a larger number of objects, 
possibly spanning a large range of redshifts, environments, and 
mass. 

Strong gravitational lensing provides the most direct probe of 
mass in early-type galaxies: a measurement of the mass enclosed 
within the Einstein radius. When combined with a spatially re- 
solved measurement of the line-of-sight velocity dispersion profile 
and the surface brightness of the lens galaxy, this mass measure- 
ment can be used to constrain the luminous and dark matter mass 
profiles simultaneously through the Jeans equation (e.g., KT). Un- 
fortunately, only a handful of known strong lenses are amenable 
to this type of analysis, often because the background source is a 
bright QSO whose images outshine the lens galaxy. 

We have therefore initiated the Sloan Lens ACS (SLACS) 
Survey (see Bolton et al. 2005) to discover a much larger sample 
of new early-type strong gravitational lenses suitable for detailed 
photometric and dynamical study and thereby realize the full po- 
tential of gravitational lensing as a probe of early-type galaxy 
structure. SLACS uses the method described by Bolton et al. 
(2004, hereafter B04) to select candidate galaxy-scale gravita- 
tional lens systems from within the Sloan Digital Sky Survey 
(SDSS) spectroscopic database on the basis of multiple higher 
redshift emission lines in the spectrum of lower redshift target 
galaxies. For a discussion of other spectroscopic lens surveys and 
lens discoveries based on the identification of anomalous emis- 
sion lines, see Huchra et al. (1985), Warren et al. (1996, 1998, 
1999), Willis (2000), Hall et al. (2000), Hewett et al. (2000), 
Johnston et al. (2003), and Willis et al. (2005). The SDSS data also 
provide stellar velocity dispersion measurements for all SLACS 
lens candidates. The most promising candidates (see § 2.1) are 
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observed through F435W and F814W using the Advanced Camera 
for Surveys (ACS) aboard the Hubble Space Telescope (HST). 
ACS images enable us to confirm bona fide lenses, measure de- 
tailed photometric and morphological parameters of the lensing 
galaxies, and obtain accurate astrometry and surface brightness 
measurements of strongly lensed features with which to constrain 
gravitational lens models. Finally, confirmed lenses from the pro- 
gram are being targeted for deep ground-based spectroscopy in 
order to obtain spatially resolved velocity dispersion profiles of 
the lensing galaxies, which will provide more precise dynamical 
constraints on their mass structure. 

This paper (Paper I) is the first in a series presenting the results 
of the current Cycle 13 HST imaging component of the SLACS 
Survey. This paper presents the catalog of new lenses confirmed 
by the SLACS HSTACS Snapshot Survey and their properties as 
measured by SDSS, describes our methods of image processing 
and lens-galaxy image subtraction, and provides an analysis of 
the possible selection biases that bear on whether or not our lenses 
are a representative sample of similar SDSS early-type galaxies. 
We also present ground-based integral-field spectroscopy of sev- 
eral SLACS systems, which supports a strong-lensing interpre- 
tation of the observed features. Treu et al. 2006 (hereafter Paper II) 
presents photometric and morphological measurements of the 
SLACS lens sample from HST imaging and places the sample within 
the context of the FP, Koopmans et al. 2006 (hereafter Paper IIT) 
presents gravitational lens and dynamical modeling results, and 
Paper IV focuses on the properties of the lensed source galaxies. 

The outline of this paper is as follows. In § 2, the survey selec- 
tion procedure and the available data of each system are described. 
In § 3, 19 newly discovered lens systems and nine unconfirmed 
systems or nondetections are presented. In § 4, the statistics and 
possible selection effects of the sample are discussed; § 5 discusses 
some of the implications of the SLACS survey for future sur- 
veys, and in § 6 conclusions are drawn and future work is dis- 
cussed. Throughout this paper, we assume a cosmological model 
with Qn = 0.3, Qa = 0.7, and Ho = 70 h7 km s7! Mpc"! (with 
hyo = 1). 


2. THE SURVEY 
2.1. Candidate Selection 


The method by which we select our lens candidates is essen- 
tially that described by B04. Briefly, we subtract best-fit principal- 
component analysis (PCA) templates (a by-product of the redshift 
pipeline) from the observed SDSS target galaxy spectra. We re- 
quire the SDSS continuum to be well fit by the template, which 
effectively yields a parent sample of galaxies with well-behaved 
absorption-dominated spectra and very secure redshifts, which 
we denote zpg (for “foreground’’). The residual spectra are then 
scanned for nebular line emission at redshifts greater than Zpg. 
Spectra in which such emission is significantly detected for at least 
three separate common atomic transitions at a single background 
redshift (zgq) are taken as lens candidates for having emission 
at two different redshifts along the same line of sight as sampled 
by the 3” diameter SDSS spectroscopic fiber. 

Taking Zp, Zgo, and the measured 3” stellar velocity dispersion 
Ga from SDSS allows us to determine to first order the strong- 
lensing cross section for each system in the source plane, using a 
singular isothermal sphere (SIS) model, as 707, with the Einstein 
radius given by 0g = An(a2/c7)(Dis/Ds). In the image plane the 
region of multiple imaging is 4 times larger than the source 
plane strong-lensing cross section, in the sense that all images 
within a radius 26, of the potential center show strongly lensed 


No. 2, 2006 SLACS. I. 705 
TABLE 1 
PROPERTIES OF OBSERVED SYSTEMS 
Regt Oa 

System Plate-MJD-Fiber gr, i (arcsec) ZEG ZBG (km s~!) Sample Lens 
SDSS J003753.21—094220.1........ 0655-52 162-392 18.00, 16.81, 16.39 2.16 + 0.06 0.1954 0.6322 265 + 10 LRG Yes 
SDSS J021652.54—081345.3........ 0668-52 162-428 19.07, 17.46, 16.90 3.05 + 0.13 0.3317 0.5235 332 + 23 LRG Yes 
SDSS J0O73728.45+321618.5......... 0541-51959-145 19.38, 17.84, 17.15 2.16 + 0.13 0.3223 0.5812 310 + 15 LRG Yes 
SDSS J081931.92+453444.8......... 0441-51868-108 18.63, 17.51, 17.07 2.32 + 0.13 0.1943 0.4462 231 + 16 MAIN ? 
SDSS J091205.30+002901.1......... 0472-51955-429 17.31, 16.22, 15.78 3.36 + 0.05 0.1642 0.3239 313 + 12 LRG Yes 
SDSS J095320.42+520543.7......... 0902-52409-577 18.57, 17.61, 17.22 1.77 + 0.09 0.1310 0.4670 207 + 14 MAIN ? 
SDSS J095629.77+510006.6......... 0902-52409-068 18.41, 17.17, 16.62 2.33 + 0.09 0.2405 0.4700 299 + 16 LRG Yes 
SDSS J095944.07+041017.0......... 0572-52289-495 18.52, 17.48, 17.02 1.21 + 0.04 0.1260 0.5350 212 + 12 MAIN Yes 
SDSS J102551.31—003517.4 0272-51941-151 17.07, 16.03, 15.57 4.05 + 0.08 0.1589 0.2764 247 + 11 LRG ? 
SDSS J111739.60+053413.9. 0835-52326-571 18.72, 17.56, 17.12 2.49 + 0.11 0.2285 0.8232 279 + 21 MAIN ? 
SDSS J120540.43+491029.3. 0969-52442-134 18.40, 17.22, 16.65 2.30 + 0.10 0.2150 0.4808 235 + 10 MAIN Yes 
SDSS J125028.25+052349.0. 0847-52426-549 18.40, 17.26, 16.77 1.76 + 0.07 0.2318 0.7946 254 + 14 MAIN Yes 
SDSS J125135.70—020805.1 0337-51997-480 18.58, 17.59, 17.24 3.64 + 0.19 0.2243 0.7843 216 + 23 MAIN Yes 
SDSS J125919.05+613408.6. 0783-52325-279 18.80, 17.46, 17.01 1.94 + 0.07 0.2333 0.4488 263 + 17 LRG ? 
SDSS J133045.53—014841.6 ; 0910-52377-503 18.34, 17.45, 17.05 0.84 + 0.04 0.0808 0.7115 178 + 09 MAIN Yes 
SDSS J140228.21+632133.5......... 0605-52353-503 18.26, 16.98, 16.49 2.67 + 0.08 0.2046 0.4814 275 + 15 LRG Yes 
SDSS J142015.85+601914.8......... 0788-52338-605 16.39, 15.56, 15.17 2.17 + 0.03 0.0629 0.5350 194 + 05 MAIN Yes 
SDSS J154731.22+572000.0......... 0617-52072-561 17.94, 16.84, 16.39 2.56 + 0.06 0.1883 0.3955 243 + 11 LRG ? 
SDSS J161843.10+435327.4......... 0815-52374-337 18.78, 17.60, 17.09 1.34 + 0.05 0.1989 0.6656 257 + 25 MAIN Yes 
SDSS J162746.44—005357.5........ 0364-52000-084 18.54, 17.29, 16.89 2.08 + 0.08 0.2076 0.5241 275 + 12 LRG Yes 
SDSS J163028.15+452036.2......... 0626-52057-518 18.84, 17.41, 16.92 2.02 + 0.07 0.2479 0.7933 260 + 16 LRG Yes 
SDSS J163602.61+470729.5......... 0627-52 144-464 18.92, 17.68, 17.18 1.48 + 0.05 0.2282 0.6745 221 + 15 MAIN ? 
SDSS J170216.76+332044.7......... 0973-52426-464 18.01, 16.91, 16.40 2.80 + 0.07 0.1784 0.4357 239 + 14 LRG ? 
SDSS J171837.39+642452.2......... 0352-51789-563 16.86, 15.97, 15.54 3.67 + 0.07 0.0899 0.7367 270 + 16 MAIN Yes 
SDSS J230053.14+002237.9......... 0677-52606-520 18.97, 17.63, 17.14 1.76 + 0.10 0.2285 0.4635 283 + 18 LRG Yes 
SDSS J230321.72+142217.9......... 0743-52262-304 17.58, 16.39, 15.96 3.02 + 0.09 0.1553 0.5170 260 + 15 LRG Yes 
SDSS J232120.93—093910.2........ 0645-52203-517 16.07, 15.21, 14.82 3.92 + 0.05 0.0819 0.5324 236 + 07 MAIN Yes 
SDSS J234728.08—000521.2........ 0684-52523-311 19.81, 18.51, 17.95 1.78 + 0.28 0.4168 0.7145 330 + 50 LRG 2 


Nortes.—The column headed “System” gives truncated J2000.0 right ascension and declination in the format HHMMSS.ss4 


tDDMMSS:s. De Vaucouleurs model 
SDSS (AB) magnitudes have been dereddened using dust maps from Schlegel et al. (1998) and have statistical errors of approximately 0.01. De Vaucouleurs 
effective radii are quoted at the intermediate axis. The column headed “Sample” indicates whether the candidate was selected from the SDSS LRG sample or from 


the SDSS MAIN galaxy sample (see § 2.1). 


counterimages. To maximize the number of strong lenses in our 
survey, we formed our Cycle 13 HST target list from the candi- 
dates with the highest predicted lensing cross section. Thus, our 
HST target sample is velocity dispersion selected to leading order, 
with an additional selection bias in favor of systems with larger 
angular diameter distance ratios D),,/D, (a function of the red- 
shifts zgg and zgc). We take 20 of our targets from the list of lens 
candidates published in B04, which were selected from within the 
luminous red galaxy (LRG) spectroscopic sample of the SDSS 
(Eisenstein et al. 2001). The LRG sample is defined by photo- 
metric selection cuts that very efficiently select massive early-type 
galaxies in the redshift range 0.15 < zrg S 0.5, as confirmed by 
SDSS spectroscopy. These galaxies are very homogeneous in their 
spectral, photometric, and morphological properties, and we place 
no further requirements on our LRG sample lens candidates be- 
yond their photometric selection as described in Eisenstein et al. 
(2001) and subsequent spectroscopic confirmation as galaxies 
with redshifts z > 0.15. The remaining 29 targets on our Snap- 
shot list are selected with the same spectroscopic algorithm from 
within the MAIN galaxy sample of the SDSS (Strauss et al. 2002). 
The MAIN sample is much more heterogeneous, and we impose 
a quiescent, absorption-dominated spectral criterion by requiring 
our lens candidates to have rest-frame equivalent widths in Ha 
of EWpa < 1.5 A (with a few exceptions made in the interest of 
maximizing the total number of discovered lenses in the program; 
see SDSS J1251—0208 below). Although most of the systems we 
target do indeed exhibit early-type morphology, the selection of 


SLACS lens candidates from within the LRG and MAIN galaxy 
samples of the SDSS is purely spectroscopic. 


2.2. ACS Image Processing 


Our ACS Snapshot observing strategy is discussed in Bolton 
et al. (2005) and consists of one 420 s Wide-Field Channel (WFC) 
exposure through each of the two filters F435W and F814W. The 
SLACS lens candidates are selected to yield bright lenses with 
faint background sources. This facilitates the photometric, mor- 
phological, and kinematic study of the lens galaxy, but can also 
make the relatively faint lensed features difficult to detect and even 
more difficult to use in lens modeling. The key to success is ef- 
fective subtraction of the image of the lens galaxy by fitting a 
smooth model to the image data of the candidate lens galaxy (e.g., 
Peng et al. 2002; Simard 1998). Given the extremely regular iso- 
photal structure of most early-type galaxies, this is a reasonable 
proposition. The most common parameterized model for early- 
type brightness distributions is the generalized de Vaucouleurs 
or Sérsic law (Sersic 1968; Ciotti 1991). We first attempted to use 
elliptical Sérsic model-fitting and subtraction in our analysis, but 
generally found that the global systematic residuals of the fit were 
large compared to the surface brightness levels of the lensed fea- 
tures that we hope to use to constrain gravitational lens models. 
Significant galaxy-core residuals are especially pronounced in the 
Sérsic-subtracted F814W data. The shortcomings of the Sérsic 
model led us to implement a more generalized galaxy model- 
fitting procedure involving a b-spline fit to the radial profile with a 
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Fic. 1.—SDSS spectroscopy showing background line emission in SLACS target galaxies. Upper gray line shows SDSS template fitted to the continuum of the 


foreground galaxy, and lower gray lines show the 1 ø noise level. 


low-order multipole dependence to fit the angular structure, which 
we describe in Appendix A. 


3. NEW LENSES AND OTHER OBSERVED SYSTEMS 


Here we present the catalog of 28 candidate lens systems? ob- 
served by the SLACS survey through 2005 March 22. Photo- 
metric and spectroscopic parameters for the sample as measured 
by SDSS are given in Table 1. Figure | presents the SDSS dis- 
covery spectra of all targeted systems, focused on the background 
redshift line emission. For definiteness, we enumerate the sev- 
eral possible explanations for any one of our spectroscopically 
selected lens candidates as follows: (1) a multiply imaged back- 


8 As of 2006 January, another 11 targets have been observed under the 
Cycle 13 SLACS program. 


ground galaxy (“strong lens”), (2) a singly imaged, although 
possibly magnified, background galaxy (a “nonlens” in our clas- 
sification), (3) a multiple-image system corresponding to a pro- 
jection of multiple, singly imaged sources with similar colors, 
and (4) spurious noise features in the SDSS spectrum. Given the 
adopted significance threshold for selection and the careful treat- 
ment of the noise properties of SDSS spectra (B04), we expect 
(4) to be an unlikely explanation, although it merits consider- 
ation in the case of any lower signal-to-noise ratio candidates 
that do not produce detections in follow-up observations. 
Systems falling into category (3) could in principle be mistaken 
for bona fide strong lenses, but can be expected to have ill- 
fitting lens models, since real lenses occupy only a small sub- 
space of all conceivable image configurations. For simplicity, 
this paper only categorizes ACS-observed systems into strong 
lenses on the one hand and everything else on the other, where 
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everything else includes systems for which a definitive expla- 
nation cannot be made based on the data in hand (i.e., ACS 
nondetections). 

To determine the incidence of strong lensing, we examine 
the F814W and F435W residual images formed by subtracting 
smooth models of the foreground galaxy constructed as described 
in Appendix A. If the residual image of a target shows multiple 
images with similar colors that can be reproduced by a simple 
lens model with the potential center fixed at the optical center of 
the foreground galaxy, we classify the system as a strong gravi- 
tational lens. The lens models for these systems are presented in 
Paper III, but we summarize the lens-modeling process here. The 
procedure is based on the method of Warren & Dye (2003), with 
an implementation as described in Treu & Koopmans (2004) and 
Koopmans (2005). We parameterize the lens galaxy as a singu- 
lar isothermal ellipsoid (SIE; Kormann et al. 1994), with three 
free parameters corresponding to lens strength (i.e., Einstein ra- 


dius), ellipticity, and major-axis position angle. The source plane 
brightness distribution of the candidate lensed galaxy is repre- 
sented on a pixel grid. We determine the best-fit lens parameters 
and source plane pixel brightnesses by minimizing x? with re- 
spect to the galaxy-subtracted /-band image-plane data, while also 
imposing a brightness-dependent regularization on the source- 
plane distribution to suppress curvature. Unsuccessful strong-lens 
models are characterized by an inability to map putative multiply 
imaged features to a single region of the source plane for any set of 
SIE model parameters. SDSS J1025 is an example of a system 
with an apparent arc and counterimage that nonetheless cannot be 
reproduced by an SIE lens model and is thus not classified as a 
strong lens. ACS data for target galaxies classified as lenses are 
shown in Figure 2. We note again that all target systems have line 
emission in their SDSS spectra at a redshift higher than that of the 
foreground galaxy, and thus, the lensing interpretation does not 
rest on imaging data alone. 
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Figure 3 presents all systems other than those that we classify 
solidly as lenses. This figure encompasses systems that may in- 
deed be strong lenses, as well as systems that appear to be definite 
nonlenses. We note that a significant number of systems that we 
consider possible but not definite lenses show a faint candidate 
counterimage near the center of the foreground galaxy opposite a 
more prominent image at a larger radius. If these features were all 
due to simple residual error in the foreground galaxy subtraction, 
this configuration would not be expected. The features could 
conceivably be explained by the prominent image “pulling” the 
foreground galaxy model to the side, leading to an undersubtrac- 
tion of the foreground galaxy flux near the center on the opposite 
side. We believe this explanation is unlikely since we have taken 
care to mask neighboring images when fitting the foreground 
galaxy model (Appendix A), and we do not see any corresponding 
oversubtraction on the near side of the core. The most asymmetric 


double lenses will in general have faint counterimages at a small 
angular offset from the lens center and will be the most difficult 
lenses to confirm, although their status as lenses or nonlenses can 
have a significant impact on statistical inferences based on the lens 
sample. Integral-field spectroscopy provides the best chance to 
associate background redshift line emission with these faint, pu- 
tative counterimages near the lens center and thereby confirm or 
reject a strong-lensing hypothesis. 

Many of the lens galaxies shown in Figure 2 show such 
striking features in their residual images that their interpretation 
as strong lenses is effectively certain when considered in combi- 
nation with their SDSS spectroscopic detection and gravitational 
lens models (Paper III). For putative lenses with less dramatic 
morphology, further evidence as to their status may be obtained 
with integral-field spectroscopy, which can confirm or deny 
the spatial coincidence between the high-redshift line emission 
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and the candidate lensed features seen with HST. For several of 
our SLACS target galaxies, we have obtained spatially resolved 
spectroscopy with the integral-field units (IFUs) of the Inamori 
Magellan Areal Camera and Spectrograph (IMACS; Bigelow & 
Dressler 2003; Schmoll et al. 2004) on the 6.5 m Walter Baade 
(Magellan I) telescope at Las Campanas Observatory and the 
8 m Gemini-North Multi-Object Spectrograph (GMOS-N; Hook 
et al. 2003; Murray et al. 2003) at Mauna Kea; we present these 
data in Figure 4. In all cases, the IFU emission-line images are 
coincident with the ACS residual-image morphology, support- 
ing the strong-lensing interpretation. 

We note here that the background sources in our lens systems 
are generally faint galaxies with irregular morphology. As evi- 
denced by the emission-line flux by which they are originally de- 
tected, they should also be regarded as star forming. The absence 


of significant broadening of the Balmer emission lines at the 
~150 km s~! resolution of the SDSS further suggests that they 
do not host quasars. For a recent observational study of the likely 
source population of emission-line galaxies (see Drozdovsky et al. 
2005). 


4. STATISTICS AND SELECTION EFFECTS 


The statistics of strong gravitational lensing are the subject 
of a great body of literature (e.g., Turner et al. 1984; Fukugita & 
Turner 1991; Kochanek 1996; Chae et al. 2002), most of it deal- 
ing with the lensing of quasars by galaxies. The rigorous statis- 
tics of the SLACS survey, although involving similar considerations, 
would be sufficiently distinct to require their own detailed treat- 
ment: our original sample is formed from potential lensing ob- 
jects rather than from potentially lensed sources, and we must 
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consider the effects of an extended source image and finite fiber 
sampling. We defer a full-lens statistical analysis of the SLACS 
survey to a future paper. Here we address the selection effects of 
the SLACS survey in terms of the robust statistics of the distri- 
bution of lens-galaxy observables. 

As discussed above, the SLACS target selection involves an 
explicit velocity dispersion and redshift selection. This should 
not in itself present any complication in the comparison of our 
lenses to other galaxies with similar velocity dispersions and red- 
shifts. After controlling for this selection, we may ask whether our 
lenses are a representative sample of the populations from which 
they have been drawn, i.e., LRGs and quiescent MAIN sample 
galaxies from the SDSS. Ifso, their relation to early-type galaxies 
selected by other means is simply determined by the SDSS LRG 
and MAIN sample target selection function (Eisenstein et al. 
2001; Strauss et al. 2002). The question can be framed between 
two alternatives: do we select simply for spectral superposition 
caused by a chance foreground-background coincidence, or is 
there some significant bias by which galaxies with certain prop- 
erties are more likely to be selected as lens candidates based on 
their SDSS spectra? To test this possibility, we exploit the parent 
sample of ~10° SDSS galaxies (Early Data Release through Data 
Release 3). For each lens we construct a control sample of SDSS 
galaxies with the same redshift and velocity dispersion (within the 
uncertainties) — quantities for which we explicitly select—and test 
whether the lens galaxy is typical of galaxies in the control sample 
in terms of its magnitude, color, effective radius, and isophotal 
ellipticity—quantities for which we do not explicitly select. Al- 
though this is not a definitive test of the representative nature of 
our lenses (it will not uncover any “hidden variable” bias), it is 
straightforward, robust, and informative. By forming control sam- 


ples at the redshift of the lenses, we avoid the necessity of ap- 
plying evolution- and k-corrections to the broadband magnitudes. 
We may perform our comparison sensibly for galaxies with ve- 
locity dispersions well measured by SDSS. We exclude the lens 
SDSS J1251—0208 from the analysis, both because its spectral 
signal-to-noise ratio is too low for a confident velocity dispersion 
measurement and because its Ha equivalent width exceeds our 
quiescent threshold. (In fact, it is a bulge-dominated spiral galaxy 
lens.) The double-lens galaxies SDSS J1618+4353 and SDSS 
J1718+6424 are also excluded, as is the lens SDSS J1205+4910, 
which has significant flux from a neighboring galaxy within the 
SDSS spectroscopic fiber. This leaves a well-defined sample of 
15 single early-type lenses. The comparison sample for each lens 
consists of unique galaxies from the SDSS database with redshifts 
within +0.005 of the lens redshift, velocity dispersions within 
+15 km s7! (approximately the median velocity dispersion error 
of the sample) of the lens galaxy, and median signal-to-noise ratio 
per spectral pixel greater than 8. For LRG sample targets, we also 
require the comparison sample to pass the photometric LRG cuts, 
and for the MAIN sample lenses we require a rest-frame EWy. < 
1.5 A. We make a luminosity-distance correction to the broad- 
band magnitudes and an angular diameter—distance correction 
to the effective radii within the redshift slice to place all galaxies 
closer to the exact redshift of the lens. The resulting individual 
comparison samples have as few as 17 and as many as 1793 gal- 
axies, with a mean of 451 and a median of 337 control galaxies 
per lens galaxy. 

We would like to answer the question of whether our lenses 
have observables (luminosities, colors, effective radii, and el- 
lipticities) consistent with having been drawn from the distribu- 
tion of those observables seen in the control samples. Since we 
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only have one lens at each redshift, a straightforward Kolmogorov- 
Smirnov (K-S) test of the lens sample against the control samples 
is not possible: even under the null hypothesis, each individual lens 
has a different parent distribution. We may, however, put all these 
parent distributions on a common footing by using the fact that the 
K-S test is invariant under a monotonic rescaling of the variable 
under consideration. Specifically, for each individual lens’s control 
sample, there exists a transformation of the observable of interest 
that converts the control distribution into a uniform distribution 
between the minimum and maximum values. The correspond- 
ingly transformed lens observable is simply equal to the normal- 
ized rank (between 0 and 1) of the lens quantity within the 
cumulative distribution of the control sample in that quantity. 
Thus, we may perform a one-dimensional K-S test of the dis- 
tribution of lens-observable ranks against a uniform distribution 
over the interval 0—1. This test is in some sense like a rank cor- 
relation test within the K-S formalism: does ranking as a lens cor- 
relate with ranking in luminosity, color, size, or ellipticity? 


Figure 5 shows the cumulative distributions of lens rankings 
in r (magnitude), g — r (color), Re (effective radius), b/a (iso- 
photal axis ratio), and J, (r-band effective surface brightness: Je œ 
1070-4" R7?) within their control samples. From these distributions 
we compute the statistic Dgs, equal to the maximum difference 
between the cumulative rank distribution and the null hypothe- 
sis distribution. Since the parent distribution under the null hy- 
pothesis is known exactly by construction (uniform probability 
of any rank between 0 and 1) and since the total number of con- 
trol galaxies is much larger than the number of lenses, we com- 
pute the significance of Dgs for a distribution of 15 values against 
a known parent distribution (e.g., Press et al. 1992). (The statis- 
tical significance of an outlying lens is limited by a smaller control 
sample in that its rank is at least 1/N,q) and at most 1 —1/Nga1, where 
Ngai is the number of galaxies in its control sample.) The resulting 
probabilities of the lens sample having been drawn at random from 
the control-sample distributions in r, g — r, Re, b/a, and I, are 
0.097, 0.655, 0.264, 0.550, and 0.085, respectively. In terms of 


712 BOLTON ET AL. Vol. 638 


SDSSJ0912+0029 


ACS—WFC F814W F814W residual —. ACS—WFC F435W 


SDSSJ0956+5100 


w 


ACS-WFC F814W F814W residual ACS-WFC F435W 
SDSSJ0959+0410 


` 


ACS-WFC F814W F814W residual ACS—WFC F435W 


SDSSJ1205+4910 


+ 


> 
7 


ACS—WFC F814W F814W residual ACS—WFC F435W F435W residual 


Fic. 2.— Continued 


color and axis ratio, the SLACS lenses seem to be a representative This trend is in the same sense as that of the FP: at fixed ve- 
sample of the parent distribution of SDSS galaxies. However, the locity dispersion, brighter galaxies have smaller effective radii. 
distribution of SLACS galaxies in brightness and effective radius Paper II shows that the SLACS lenses fall on the local FP cor- 
(and hence in surface brightness) is somewhat significantly skewed rected for stellar evolution and further examines the degree of lens 


toward brighter and more centrally concentrated systems. bias toward higher effective surface brightness in the context of 


No. 2, 2006 


SDSSJ1250+0523 


ACS—WFC F814W 


SDSSJ1251—-0208 


ACS—WFC F814W F814W residual 


SDSSJ 1330-0148 


ACS-WFC F814W 


SDSSJ1402+6321 


ACS-WFC F814W F8t4W residual- 


SLACS. I. 713 


ACS—WFC -F435W 


` 


ACS-WFC F435W 


F435W- residual 


ACS-WFC F435W 


ACS-WFC F435W 


Fic. 2.— Continued 


the FP. Although the statistical significance is not great, the trend 
is at first glance suggestive of the effect that more centrally con- 
densed objects are more efficient gravitational lenses (e.g., Li 
& Ostriker 2002), to the extent that stellar mass constitutes the 
dominant lensing component. A larger sample of similarly selected 
lenses will allow us to test this selection effect with greater sta- 


tistical significance. Possible interpretations of this finding—as 
discussed in detail in Papers II and II]—include observational 
selection effects due to the finite size of the SDSS fiber, which 
may bias our survey toward the highest surface brightness lens 
galaxies. In Paper III we will see that, within the context of 
power-law mass models, the SLACS lens sample shows great 
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homogeneity in the slope of the total mass density profile (y’) 5. DISCUSSION AND IMPLICATIONS 

and no significant correlation between y’ and J,. This uniform FOR FUTURE SURVEYS 

degree of mass concentration within the lens sample itself argues 

somewhat against a lens selection bias due to mass concentration, The large number of new strong lenses confirmed by the 
but this argument is limited by the lack of an equally direct probe SLACS Survey demonstrates the technical feasibility of carry- 


of the mass distribution in the nonlens control sample. ing out spectroscopic surveys for strong gravitational lenses and 
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suggests that similar gravitational lens surveys should be consid- 
ered as explicit goals of future large-scale spectroscopic surveys. 
The SLACS results are particularly notable when one considers 
that this type of lens survey was never one of the science or design 
goals of the SDSS, and the survey is thus a testament to the 
performance of the SDSS telescope and spectrograph and to the 


optimal quality of the SDSS spectroscopic data pipeline. With the 
clarity of hindsight, we may identify those particular features of 
the SDSS that were beneficial to the spectroscopic lens survey 
project and those that could conceivably be improved in any 
future survey to increase the yield of spectroscopically selected 
strong lenses suitable for specified science goals. 
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Fic. 3—HST ACS-WFC imaging of lens candidates observed by the SLACS survey for which the ACS data are inconclusive as to the incidence of strong lensing. 


Images are as in Fig. 2. See notes on individual systems in Appendix B. 


The large parent sample of galaxies observed is a clear ad- 
vantage afforded by the SDSS for this project or any other rare 
object survey. Empirically, we have found that spectroscopic 
lenses and lens candidates appear with a frequency on the order of 
1 in a 1000 in the SDSS; starting with ~105 early-type galaxies is 
therefore one key to success. At the same time, the large number 
of highly magnified galaxies that we have found close to the 
spectroscopic detection limit (note the generally low emission-line 
signal-to-noise ratio in the discovery spectra of Fig. 1) suggests 
that we may only be seeing the tip of the iceberg. We conclude that 
(1) a spectroscopic lens survey will always be at an advantage if it 
can piggyback on a more broadly conceived spectroscopic survey 
that does not have gravitational lens discovery as its sole scientific 
goal and (2) future spectroscopic surveys that go significantly 
deeper than the SDSS should, all other factors being equal, dis- 
cover a significantly larger fraction of gravitational lenses. 

One might suspect that a significant factor for the SLACS sur- 
vey is the large (3” diameter) aperture of the SDSS spectroscopic 
fibers. Certainly, a larger fiber will always collect more photons, 


but for lenses with image splittings of ~2” (i.e., ~1” Einstein 
radii), such as we have presented here, a 3” fiber could afford an 
advantage over smaller fibers disproportionate to its larger aper- 
ture due to the concentration of lensed flux away from the center 
of the fiber. We can test the importance of this aperture size ef- 
fect by simulating observations with varying fiber apertures. We 
assume an Einstein radius of 1”2 (an approximate median value 
from Paper III) and an exponential disk background galaxy with a 
disk scale length of 0”2. We compute lensed images for a range of 
source galaxy impact parameters using an SIS lens model. We then 
smear the images with 173 seeing (the median SDSS spectro- 
scopic seeing for our targets) and integrate over the 3” diameter 
SDSS fiber. If we reduce the fiber aperture to 2” in diameter, we 
find that the flux collected from the lensed galaxy is reduced to 
a fraction of 0.4—0.5, depending on impact parameter. This is 
essentially the factor of (1/1.5)? ~ 0.44 corresponding to the di- 
minished aperture size. Thus, for lenses of the angular size typical 
of SLACS lenses, the advantage of a large fiber is just the simple 
advantage of a large aperture due to the smearing effects of seeing; 
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a smaller fiber would not have lost a disproportionate amount of 
line flux due to lensing geometry. However, for lenses with larger 
image splittings (which would occur for higher redshift sources 
behind the same foreground galaxies), a small fiber would likely 
be a qualitative disadvantage to a spectroscopic lens survey. Be- 
cause a fiber diameter much greater than 3” would probably in- 
crease spurious nonlensed interlopers, we conclude that the SDSS 
fiber size represents a very good compromise when searching for 
massive early-type lens galaxies at these redshifts. 

For two principal reasons, our spectroscopic survey owes 
much success to having targeted lens candidates with background 
redshifts confirmed by multiple emission lines. First, the inci- 
dence of false-positive spectroscopic detections and emission- 
line misidentifications is negligible. Second, to fixed line flux 
limits, intermediate-redshift emission from [O un] 3727, [O n] 
5007, and the hydrogen Balmer series is more common than Lya 
emission at high redshift (Hippelein et al. 2003; Maier et al. 
2003), and thus, we see an abundance of oxygen/Balmer lenses. 
However, for optical surveys, this survey strategy limits the back- 
ground redshifts to zgg < 0.8, beyond which HG and [O m] 5007 


move out of the observable band. There is a definite incentive 
to discover significant numbers of lenses with higher lens and 
source redshifts, both to probe evolution in the lens population 
and to observe lensed images at larger physical radii within the 
lens galaxy in order to obtain greater leverage on dark matter halos. 
Therefore, in designing future surveys to discover gravitation- 
ally lensed emission-line galaxies at higher redshift, two obvious 
considerations are increased survey depth to detect a fainter popu- 
lation and increased spectroscopic resolution beyond the 1/Ad = 
1800 resolution of the SDSS in order to split the [O n] 3727 
doublet and resolve the characteristic asymmetric profile of Lya 
emission (thus permitting more secure single-line redshifts). 
One of the greatest uncertainties at the outset of the SLACS 
survey was simply the fraction of strong lenses within the can- 
didate list. One can attempt to calculate the lensing probabilities 
of the candidates given Zpg, zga, and og, assuming a model for 
the SDSS spectroscopic observations and assuming a luminosity 
function (LF) and size/shape for the background emission-line 
galaxies, as was done in B04, but the ingredients are all uncer- 
tain. To put it simply, we did not know beforehand whether or 
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not the majority of our targets would owe their high-redshift line 
emission to the point-spread function (PSF)—smeared wings of 
bright, unlensed galaxies at large impact parameter rather than 
to Einstein ring images of faint galaxies closely aligned with the 
SDSS target galaxy. The high ratios of lenses to nonlenses and 
of quads/rings to doubles that we see are suggestive of a large 
magnification bias in our selection, with highly magnified faint, 
lensed galaxies detected with greater frequency than less mag- 
nified lensed galaxies or unlensed projections. Future work will 
quantify the extent of any such magnification bias and will derive 
statistical constraints on the lens and source populations based 
on the incidence of lensing within the survey. 


6. CONCLUSIONS AND FUTURE WORK 


The HST Cycle 13 SLACS Survey’ has produced a catalog of 
19 previously unknown early-type strong gravitational lens gal- 
axies. The >68% fraction of genuine lenses within the first 28 can- 
didates observed shows that the survey strategy—spectroscopic 


° HST program 10174; PI: Koopmans. 


candidate selection from within the SDSS, followed by HST 
ACS Snapshot observations of the systems with the largest es- 
timated lensing cross sections—is an efficient means to discover 
new gravitational lenses; similar lens surveys should be consid- 
ered an explicit goal of future spectroscopic surveys. By tar- 
geting galaxies with the highest estimated lensing cross sections, 
the Cycle 13 SLACS survey has effectively selected high-mass 
lensing galaxies (strong lensing cross section is simply propor- 
tional to lens mass within the Einstein radius). The natural selec- 
tion of lenses in other surveys is also weighted by lensing cross 
section, and thus, strong lenses are typically massive galaxies. 
During HST Cycle 14, we will extend our survey to galaxies of 
lower mass in order to use strong lensing and stellar dynamics to 
measure the mass dependence of structure and mass-to-light 
ratio in early-type galaxies.'° These galaxies will have a lower 
lensing cross section and hence should have a lower lensing rate, 
although magnification bias may skew the distribution of can- 
didates in favor of lenses over nonlenses. The final combined 


10 HST program 10587; PI: Bolton; 118 Snapshot targets. 
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Fic. 4.—IFU narrowband imaging of several SLACS target systems. Left, HST-ACS F814W residual images, for reference; center, IFU emission-line images; right, 
IFU continuum images. IFU images are constructed from IFU spectra by fitting linear continuum and Gaussian emission-line models as described in A. S. Bolton & 
S. Burles (2006, in preparation). Data are from IMACS (SDSS J0037, SDSS J0216, and SDSS J2321) and GMOS-N (others). Emission-line images for SDSS J0216, SDSS 
J0956, and SDSS J1547 have been smoothed spatially with a seven-lenslet hexagonal kernel to suppress noise. Narrowband images are formed at the redshifted wave- 
length of the following background galaxy emission lines: [O m] 45007 (420037, 0737, 1402, and 1702), H8 (20956), and [O n] 243727 (A20216, 1547, 1630, and 2321). 


lens sample will be a unique resource for the detailed mea- 
surement of the mass profile of early-type galaxies within the 
effective radius. 

The SLACS sample represents the largest single catalog of 
uniformly selected early-type gravitational lens galaxies assem- 
bled to date. The lensed images provide an aperture mass con- 
straint within a typical scale of ~Rg/2. The SLACS lens galaxies 
are all brighter than their lensed background galaxies by many 


magnitudes and are thus ideally suited to detailed photometric 
measurement as with the current ACS data. They are also excel- 
lent candidates for spatially resolved spectroscopy to constrain dy- 
namical galaxy models in combination with the mass constraints 
from strong lensing. In fact, nearly all of the SLACS galaxies 
already have well-measured luminosity-weighted velocity dis- 
persions inside the 3” SDSS fiber (see Table 1), and thus, we have 
already nearly tripled the number of known gravitational lenses 
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with known stellar velocity dispersions. Paper II will present the 
constraints on the logarithmic slope of the total radial density 
profile (luminous plus dark matter) that can be obtained by com- 
bining lens models with the aperture-integrated SDSS velocity 
dispersions. We may derive tighter constraints and a significant 
decomposition into luminous and dark components by obtaining 
spatially resolved line-of-sight velocity dispersion measurements 
within the lens galaxies (see KT). 

An analysis of the distribution of SLACS lens-galaxy photo- 
metric and structural parameters within control samples from 


the SDSS database shows that the SLACS lenses are typical 
of their SDSS parent sample with regard to color and ellipticity. 
The location of SLACS lens galaxies within the FP is the sub- 
ject of Paper II. However, the SLACS lens galaxies exhibit a 
somewhat significant bias toward brighter and higher surface 
brightness galaxies at fixed redshift and velocity dispersion 
(§ 4; see also Paper II). The significance of this bias can be 
confirmed with a larger sample of similarly selected lenses, which 
we plan to obtain with the continuation of our ACS survey dur- 
ing HST Cycle 14. 
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Results in this paper are based in part on data from the Sloan 
Digital Sky Survey (SDSS) archive. Funding for the creation 
and distribution of the SDSS Archive has been provided by the 
Alfred P. Sloan Foundation, the Participating Institutions, the 
National Aeronautics and Space Administration, the National 
Science Foundation, the US Department of Energy, the Japanese 
Monbukagakusho, and the Max Planck Society. The SDSS Web 
site is http://www.sdss.org. The SDSS is managed by the Astro- 
physical Research Consortium (ARC) for the Participating In- 
stitutions. The Participating Institutions are the University of 
Chicago, Fermilab, the Institute for Advanced Study, the Japan 


Participation Group, Johns Hopkins University, Los Alamos 
National Laboratory, the Max-Planck-Institute for Astronomy 
(MPIA), the Max-Planck-Institute for Astrophysics (MPA), New 
Mexico State University, the University of Pittsburgh, Princeton 
University, the United States Naval Observatory, and the Uni- 
versity of Washington. These results are also based in part on ob- 
servations obtained with the 6.5 m Walter Baade telescope of the 
Magellan Consortium at Las Campanas Observatory. These results 
are also based in part on observations obtained under program GN- 
2004A-Q-5 at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., 
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Fic. 5.—Kolmogorov-Smimov tests of the rank of lens-galaxy observables within the cumulative distributions of those observables in control-galaxy samples with 
similar redshifts and velocity dispersions. Top left to bottom right: Distributions for r-band magnitude, g — r color, effective (half-light) radius R,, isophotal axis ratio 
b/a, and effective surface brightness Z, (all as measured from SDSS imaging data for a de Vaucouleurs galaxy model). Solid black lines show the cumulative distribution 
in rank, solid gray lines show the null hypothesis for this distribution, and vertical dash-dotted lines show the location of the maximum difference between the two 
distributions. 
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APPENDIX A 


b-SPLINE GALAXY MODEL SUBTRACTION 


The b-spline technique is a well-known method for fitting a piecewise-defined polynomial of arbitrary order to describe the 
dependence ofa series of data values on an independent variable (e.g., de Boor 1977). The coefficients of the polynomial change at break 
points in the independent-variable domain, whose spacing may be chosen to allow more or less freedom depending on the level of detail 
to be fit. A b-spline of order n (where by convention n = 1 is piecewise constant, n = 2 is linear, n = 3 is quadratic, n = 4 is cubic, and 
so on) has continuous derivatives to order n — 2 across the break points. For our galaxy modeling we use b-splines of order n = 4, which 
have continuous zeroth, first, and second derivatives. The coefficients of the b-spline are determined by fitting to the data in a least- 
squares sense, given the break point spacing and derivative continuity conditions. A b-spline model may be reinterpreted in terms of a 
number of localized basis functions within the domain, whose shapes are set by the order of the b-spline and by the break point spacing 
and whose amplitudes are determined by the fit to the data. The fit itself is entirely linear, and the localized influence of each basis 
function within the domain implies that only a banded-diagonal matrix need be inverted in the solution for the coefficients. 
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The radial b-spline technique for galaxy images permits smooth fitting of arbitrary radial brightness profiles. Since we do not have 
multiple dithered exposures to combine, we perform our b-spline galaxy model fitting in the native ACS pixel coordinates of the images. 
We account for the distortion in the ACS by using the solution provided in the image headers to compute relative tangent plane right 
ascension and declination values for all pixels in the image, which are taken as the independent variables for the fit. We fit for surface 
brightness as a function of position using cosmic-ray—masked flat-fielded images (i.e., calibrated by the CALACS software pipeline’! to 
measure surface brightness rather than counts per pixel), so our fits are not biased by slight variations in pixel area across the images. The 
one nonlinear step in our b-spline galaxy model fitting is the determination of the center of the lensing galaxy; we describe our centering 
method further below. The model-fitting procedure is carried out using the adopted center for the lens galaxy within a suitable subsection 
of the ACS field (typically 12” x 12”), along with an error image, a mask specifying cosmic-ray and other zero-weight pixels, and a 
mask corresponding to stars, neighboring galaxies, and any apparent background galaxy features. The cosmic-ray masks are initially 
generated by the LACOSMIC software (van Dokkum 2001) and adjusted manually over the image subsection surrounding the lens; the 
neighboring-object masks are created manually. For each pixel in the subimage, the radial offset R from the galaxy center (in arcseconds) 
and the azimuthal angle 6 relative to a fixed position angle are computed. A set of break points in R is chosen for the fit, typically every 
02-073 in the central 1”—2”, with increased spacing farther out. (We note again that the ACS pixel scale is 0705 per pixel.) For the case 
of complete Einstein rings, break points at the ring radius can be removed to allow smoother interpolation of the lens-galaxy model over 
the ring region. A multipole angular dependence is incorporated into the fit as follows. If the purely one-dimensional b-spline fit is 
represented as 


IR) =X ax fel), (Al) 
k 


where the f,(R) are the localized basis functions that are only nonzero over a small range of break points and the a; are their amplitudes, 
then the two-dimensional fit is represented by 


ICR, 0) = X [Bmx cos (mB) + cme sin (mB)] fe(R). (A2) 
m,k 


The number of multipole orders to be fit is chosen individually for each galaxy, with m = 0 (monopole) and m = 2 (quadrupole) 
always present and higher order terms added if necessitated by systematic angular structure in the residual images. With the occasional 
exception of an m = | (dipole) term to capture slightly disturbed morphologies or mild miscentering in the central regions, only even 
multipole orders are used. This form of angular dependence is ideal for fitting early-type galaxies: the fit remains linear, the global 
symmetries seen in early-type galaxies are naturally reflected in the low-order terms, and effects such as isophotal twists, varying 
ellipticity with radius, and diskiness/boxyness can be captured with minimal effort (with only the last of these effects requiring 
multipole orders beyond the quadrupole). We note also that the monopole term can automatically include a fit to the sky background. 

Following the completion of the initial b-spline model-fitting step (which for our 12” x 12”, 240 x 240 pixel images generally takes 
2-3 s ona 2.53 GHz Pentium 4 Linux PC), we examine the residual (data minus model) image for faint features not associated with the 
smooth galaxy model and not initially masked and perform a second fit with an updated mask, also adding multipole orders as needed. 
Since at this stage we are principally concerned with generating high-quality residual images and not with the measurement of lens- 
galaxy parameters, we fit directly to the images without convolving the b-spline model with the HST PSF. As an example, we describe 
the fit to the sky-subtracted image of the E3 lens galaxy SDSS J0912+0029 to give the F814 residual image shown in Figure 2, which 
uses break points spaced every 073 in R and multipole terms of order 0, 2, and 4 (monopole, quadrupole, and octopole). The ratio of 
the m = 2 amplitudes [defined as (b3, + c3,)!*] to the m = 0 amplitudes rises from ~0.25 in the center to ~0.4 at the effective radius 
(about 3”), and the ratio of the m = 4 to the m = 0 amplitudes rises from ~0 to ~0.1 over the same range. Some of this increase in 
angular structure with radius is due to our fitting directly to the PSF-blurred data. 

The centering necessary for the b-spline galaxy fit is accomplished by first fitting with respect to a best-guess center with monopole, 
dipole, and quadrupole terms. This fit is then evaluated in a ring and used to determine a flux centroid, which is adopted as a new best- 
guess center for another iteration of the same procedure. Within a few iterations this process—which amounts to minimizing the dipole 
term—converges to a position that we adopt as the constrained center for the final model fit. 


APPENDIX B 
NOTES ON INDIVIDUAL SYSTEMS 
B1. LENSES 


SDSS J0216—0813.—A faint but definite counterimage is seen to the northeast in both bands in this system, opposite the more 
obvious extended arc. 

SDSS JO912+ 0029.—This system shows a faint, extended counterarc to the south in the F814W band, opposite the more prominent 
northern arc. Both images are only marginally detected in the F435W band. 

SDSS J0956+ 5100.—This system includes a small round red companion to the north of the lens galaxy, not to be confused with the 
prominent lensed images. 

SDSS J1251—0208—This galaxy has extended spiral structure in addition to the lensing bulge. 


1l See http://www.stsci.edu/hst/acs. 


724 BOLTON ET AL. 


SDSS J1330—0148—This system is judged to be a lens on the basis of the compact counterimage detected at high significance in 


both bands. 


SDSS J1402+ 6321.—This lens, with faint quadruple images, was the first system observed by the SLACS survey and is the subject of 


Bolton et al. (2005). 


SDSS J1618+ 4353.—The angular resolution of HST reveals this system to be a pair of foreground galaxies. The pair lenses a com- 
pact background source into a 3+1 quadruple configuration. The extra compact image to the west of the counterimage is more blue in 
color than the lensed images and thus does not pertain to the lensed configuration. 

SDSS J1718+6424.—A second nearby galaxy contributes significantly to the lensing potential in this system. 


B2. OTHER SYSTEMS 


SDSS J1117+ 0534, SDSS J1259+ 6134.—The imaging depth is not great enough to unambiguously show whether the putative lensed 
features appear in both bands, but the general geometry is plausibly similar to J0912+0029, which has an extremely solid gravitational 


lens model (Paper III). 


SDSS J1636+4707.—The curvature of the arcs is suggestive of ring geometry, and the knot-counterknot features along the north- 
south axis may be images of each other. Deeper observations would reveal whether or not this system is a real Einstein ring, but it is not 


certain from the current data. 


SDSS J1702+ 3320.—Both F435W and F814W data show possible counterimage features with plausible relative orientations and 


separations at low signal-to-noise ratio. 
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